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Abstract 



In this paper, we use the type la supernovae data to constrain the holo- 
graphic dark energy model proposed by Li. We also apply a cosmic age test to 
this analysis. We consider in this paper a spatially flat Friedmann-Robertson- 
Walker Universe with matter component and holographic dark energy com- 
ponent. The fit result shows that the case c < 1 (c = 0.21) is favored, which 
implies that the holographic dark energy behaves as a quintom-type dark 
energy. Furthermore, we also perform a joint analysis of SNe+CMB+LSS to 
this model; the result is well improved, and still upholds the quintom dark en- 
ergy conclusion. The best fit results in our analysis are c = 0.81, = 0.28, 
and h = 0.65, which lead to the present equation of state of dark energy 
wo = —1.03 and the deceleration/acceleration transition redshift zt = 0.63. 
Finally, an expected SNAP simulation using ACDM as a fiducial model is 
performed on this model, and the result shows that the holographic dark en- 
ergy model takes on c < 1 (c = 0.92) even though the dark energy is indeed 
a cosmological constant. 



I. INTRODUCTION 



The type la supernova (SN la) [1,2] observations provide the first evidence for the ac- 
celerating expansion of the present Universe. To explain this accelerated expansion, the 
Universe at present time is viewed as being dominated by an exotic component with large 
negative pressure referred to as dark energy. A combined analysis of cosmological obser- 
vations, in particular, of the WMAP (Wilkinson Microwave Anisotropy Probe) experiment 
[3-5], indicates that dark energy occupies about 2/3 of the total energy of the Universe, 
and dark matter about 1/3. The most obvious theoretical candidate of dark energy is the 
cosmological constant A [6] which has the equation of state w — — 1. An alternative pro- 
posal is the dynamical dark energy [7,8] which suggests that the energy form with negative 
pressure is provided by a scalar field evolving down a proper potential. The feature of this 
class of models is that the equation of state of dark energy w evolves dynamically during 
the expansion of the Universe. However, as is well known, there are two difficulties arise 
from all these scenarios, namely the two dark energy (or cosmological constant) problems 
- the fine-tuning problem and the "cosmic coincidence" problem. The fine-tuning problem 
asks why the dark energy density today is so small compared to typical particle scales. The 
dark energy density is of order 10~ 47 GeV 4 , which appears to require the introduction of a 
new mass scale 14 or so orders of magnitude smaller than the electroweak scale. The second 
difficulty, the cosmic coincidence problem, states "Since the energy densities of dark energy 
and dark matter scale so differently during the expansion of the Universe, why are they 
nearly equal today"? To get this coincidence, it appears that their ratio must be set to a 
specific, infinitesimal value in the very early Universe. 

Recently, considerable interest has been stimulated in explaining the observed dark en- 
ergy by the holographic dark energy model. For an effective field theory in a box of size 
L, with UV cut-off A c the entropy S scales extensively, S ~ L 3 A^. However, the peculiar 
thermodynamics of black hole [9] has led Bekenstein to postulate that the maximum en- 
tropy in a box of volume L 3 behaves nonextensively, growing only as the area of the box, 
i.e. there is a so-called Bekenstein entropy bound, S < Sbh = tcM^L 2 . This nonextensive 
scaling suggests that quantum field theory breaks down in large volume. To reconcile this 
breakdown with the success of local quantum field theory in describing observed particle 
phenomenology, Cohen et al. [10] proposed a more restrictive bound - the energy bound. 
They pointed out that in quantum field theory a short distance (UV) cut-off is related to a 
long distance (IR) cut-off due to the limit set by forming a black hole. In the other words, if 
the quantum zero-point energy density px is relevant to a UV cut-off, the total energy of the 
whole system with size L should not exceed the mass of a black hole of the same size, thus 
we have L 3 px < LM 2 this means that the maximum entropy is in order of S B ^. When we 
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take the whole Universe into account, the vacuum energy related to this holographic prin- 
ciple [11] is viewed as dark energy, usually dubbed holographic dark energy. The largest IR 
cut-off L is chosen by saturating the inequality so that we get the holographic dark energy 
density 



take L as the size of the current Universe, for instance the Hubble scale if -1 , then the dark 
energy density will be close to the observed data. However, Hsu [12] pointed out that this 
yields a wrong equation of state for dark energy. Li [13] subsequently proposed that the IR 
cut-off L should be taken as the size of the future event horizon 



then the problem can be solved nicely and the holographic dark energy model can thus 
be constructed successfully. Some speculations on the deep reasons of the holographic dark 
energy were considered by several authors [14]; further studies on this model see also [15-23]. 
In addition, it is necessary to discuss about the choice L = H^ 1 . Even though this choice was 
argued to be unsuitable due to that it may lead to the holographic dark energy tracks the 
matter density, this does not mean that the formalism M 2 H 2 cannot be made compatible 
with the observation. There are other contexts in quantum field theory where one can 
have a dark energy behaving as M 2 H 2 without introducing the holographic principle. For 
instance, Refs. [24] nicely introduce the M 2 H 2 law from general arguments in quantum 
field theory. Actually, the first time where this law was introduced in quantum field theory 
was in the context of the renormalization group models of the cosmological constant; see 
e.g. [25]. Otherwise, within the holographic model framework, the choice L = H^ 1 can 
also be favored, if one introduce some interaction between dark energy and dark matter 
[23] . However, in this paper we restrict our attention to the holographic dark energy model 
proposed by Li [13]. 

In this paper, we will see what constraints to the holographic dark energy model are 
set by present and future SNe la observations. Recently, some constraints from SNe la on 
related model where obtained in Refs. [16,17]. In this paper, we extend the analysis carried 
out in Ref. [16]. The work presented here differs from [16] (and [17]) in the following aspects: 
(a) We not only constrain the model by means of the SNe la observations, but also test the 
fit results by using the age of the Universe; (b) When performing the analysis of the SNe 
data, we also test the sensitivity to the present Hubble parameter H in the fit; (c) For 
improving the fit result, we combine the current SNe la data with the cosmic microwave 



p x = 3c 2 M 2 L- 2 



(1) 



where c is a numerical constant, and M p = l/\^8nG is the reduced Planck mass. If we 




(2) 
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background (CMB) data and the large-scale structure (LSS) data to analyze the model; (d) 
We also investigate the predicted constraints on the model from future SNe la observations. 



II. THE HOLOGRAPHIC DARK ENERGY MODEL 

The holographic dark energy scenario may provide simultaneously natural solutions to 
both dark energy problems as demonstrated in Ref. [13]. In what follows we will review 
this model briefly and then constrain it by the type la supernova observations. In addition, 
we will also apply a joint analysis of SNe+CMB+LSS data to this model. Consider now a 
spatially flat FRW (Friedmann-Robertson- Walker) Universe with matter component p rn (in- 
cluding both baryon matter and cold dark matter) and holographic dark energy component 
px, the Friedmann equation reads 

3M*H 2 = p m + p x , (3) 

or equivalently, 

^ = n^ a - s + nJ£. (4) 

Note that we always assume spatial flatness throughout this paper as motivated by inflation. 
Combining the definition of the holographic dark energy (1) and the definition of the future 
event horizon (2), we derive 

r°° din a' c 



f 

J a 



Ha' HaVQx 
We notice that the Friedmann equation implies 



(5) 



Substituting (6) into (5), one obtains the following equation 



jH e*'/yi _ n xdx ' = ce°' 2 ^±- - 1 , (7) 

where x = In a. Then taking derivative with respect to x in both sides of the above relation, 
we get easily the dynamics satisfied by the dark energy, i.e. the differential equation about 
the fractional density of dark energy, 

n' x = n x (i-n x )(i + *yfite) , (8) 

where the prime denotes the derivative with respect to x. This equation describes the 
behavior of the holographic dark energy completely, and it can be solved exactly [13,16], 
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lnn x - -^—\n(l - Jr^) + _£_ln(l + J^)-_ -ln(c + 27^) = - \n(l + z) + y , 

I + c v 2 — c v 4 — c z v 

(9) 

where yo can be determined through (9) by replacing fix with fl x as z = 0. From the energy 
conservation equation of the dark energy, the equation of state of the dark energy can be 
expressed as 

w = —1 — - — . (10) 

3 d\na y ' 

Then making use of the formula p x = 1 ^n x Pm a ~ 3 anc ^ the differential equation of fl x (8), 

the equation of state for the holographic dark energy can be given [13,16,17] 

w = -hi + -^) . (ii) 

o C 

We can also give the deceleration parameter q = — a/aif 2 , in terms of fix, 

1 U V- (12) 



fix 

2 2c 



It can be seen clearly that the equation of state of the holographic dark energy evolves 
dynamically and satisfies — (l + 2/c)/3 < w < —1/3 due to < fix < 1- In this sense, 
this model should be attributed to the class of dynamical dark energy models even though 
without quintessence scalar field. The parameter c plays a significant role in this model. If 
one takes c = 1, the behavior of the holographic dark energy will be more and more like 
a cosmological constant with the expansion of the Universe, and the ultimate fate of the 
Universe will be entering the de Sitter phase in the far future. As is shown in Ref. [13], 
if one puts the parameter fl x = 0-^3 into (11), then a definite prediction of this model, 
Wq = —0.903, will be given. On the other hand, if c < 1, the holographic dark energy will 
behave like a quintom-type dark energy proposed recently in Ref. [26] , the amazing feature of 
which is that the equation of state of dark energy component w crosses the phantom divide 
line, —1, i.e. it is larger than —1 in the past while less then —1 near today. The recent fits to 
current SNe la data with parametrization of the equation of state of dark energy find that 
the quintom-type dark energy is mildly favored [27-29]. Usually the quintom dark energy 
model is realized in terms of double scalar fields, one is a normal scalar field and the other 
is a phantom-type scalar field [30,31]. However, the holographic dark energy in the case 
c < 1 provides us with a more natural realization for the quintom picture. While, if c > 1, 
the equation of state of dark energy will be always larger than —1 such that the Universe 
avoids entering the de Sitter phase and the Big Rip phase. Hence, we see explicitly, the 
determination of the value of c is a key point to the feature of the holographic dark energy 
as well as the ultimate fate of the Universe. 
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FIG. 1. An illustrative example for the holographic dark energy model. The evolutions of 
the deceleration parameter q and the equation of state of dark energy w. In this case, we take 
Q° m = 0.27. 

As an illustrative example, we plot in Fig.l the evolutions of the deceleration parameter q 
and the equation of state of dark energy w in cases of c = 1.2, 1.0, 0.8 and 0.6, respectively. 
It is easy to see that the equation of state of dark energy crosses — 1 as c < 1, which is 
in accordance with the plots of the model- independent analysis in Ref. [27]. It should be 
pointed out that the variable cosmological constant model can also give rise to a quintom 
behavior, i.e. the effective equation of state produced by this model can also cross —1 [32]. 

In the forthcoming sections, we will see what constraints to the model described above 
are set by present and future SNe la observations. In the fitting, we use the recent new high 
red-shift supernova observations from the HST/GOODS program and previous supernova 
data, and furthermore we test the fit results by means of the cosmic age data. We find 
that if we marginalize the nuisance parameter h, the fit of the SNe observation provides 
0.09 < c ^ 0.62 (la); this means the holographic dark energy behaves as a quintom. 
Further, we find that the allowed range of model parameters depends on h evidently, for 
instance, we obtain (in la) 0.21 < c < 1.17 for h = 0.64, 0.10 ^ c ^ 0.43 for h = 0.66, and 
0.03 p$ c ^$ 0.07 for h = 0.71. The allowed regions of model parameters become evidently 
smaller for larger values of h. The fit values of the matter density in this model are 
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apparently larger than the fit value appears in the ACDM (WMAP result). Now that the 
SNe data analysis evidently depends on the value of h, an important thing we should do is 
to find some observational quantities which do not depend on h to be useful complements 
of the SNe data set to probe the property of holographic dark energy. Such quantities can 
be found in CMB and LSS (1Z and A, respectively, see section IV). A combined analysis 
of SNe+CMB+LSS shows that the confidence region evidently shrinks, and the value of c 
is changed considerably, namely in la, 0.65 < c < 1.04. We will discuss the cosmological 
consequences come from the fits, and compare the case of the joint analysis with the case 
of the SNe only in detail. Interestingly, an expected SNAP simulation using ACDM as a 
fiducial model shows that the holographic dark energy model will also favor c < 1 even 
though the Universe is indeed dominated by a cosmological constant. 



III. THE CURRENT TYPE IA SUPERNOVA CONSTRAINTS 

We now perform the best fit analysis on our holographic dark energy model with data 
of the type la supernova observations. The luminosity distance of a light source is defined 
in such a way as to generalize to an expanding and curved space the inverse-square law of 
brightness valid in a static Euclidean space, 

/ C \ 1/2 r z Hz' 

Ms^) - H ^ l+ ^Lw)' (13) 

where C is the absolute luminosity which is a known value for the standard candle SNe 
la, T is the measured flux, Hq 1 (here we use the natural unit, namely the speed of light 
is defined to be 1) represents the Hubble distance with value Hq 1 = 2997.9/i~ 1 Mpc, and 
E(z) = H(z)/H can be obtained from (4), expressed as 

note that the dynamical behavior of Qx is determined by (8). The observations directly mea- 
sure the apparent magnitude m of a supernova and its red-shift z. The apparent magnitude 
m is related to the luminosity distance 4 of the supernova through 

m(z) = M + 5 log 10 (d L (^)/Mpc) + 25 , (15) 

where M is the absolute magnitude which is believed to be constant for all type la super- 
novae. The numerical parameter c of the model, the density parameter and the nuisance 
parameter h can be determined by minimizing 

^,2 _ [/fobs(^i) ~ Vth(zj)] 2 
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where the extinction-corrected distance moduli fj,(z) is defined as fj,(z) — m(z) — M, and Oi 
is the total uncertainty in the observation. The likelihood C oc e^ x ^ 2 if the measurement 
errors are Gaussian. In our analysis, we take the 157 gold data points listed in Riess et 
al. [33] which includes recent new 14 high redshift SNe (gold) data from the HST/GOODS 
program. The results of our analysis for the holographic dark energy model are displayed 
in Fig. 2. In this figure we marginalize over the nuisance parameter h and show 68%, 95% 
and 99% confidence level contours, in the (c, Jl^J-plane. The best fit values for the model 
parameters are: h = 0.66, tt° m = 0.47lg;Sf, and c = 0.21±g;^ with X 2 m - m = 173.44 (similar 
results see also [16]). We see clearly that the fit values of this model are evidently different 
from those of ACDM, i.e. the value of h is slightly smaller and evidently larger (The 
WMAP results for ACDM model are [4,34]: h = 0.71^3 and n° m = 0.27±0.04). We notice 
in this figure that the current SNe la data do not strongly constrain the parameters and 
c (in 2a), in particular c, in the considered ranges. Other observations may impose further 
constraints. For instance, the CMB and LSS data can provide us with useful complements 
to the SNe data for constraining cosmological models. But, first, we will apply a cosmic age 
test to the SNe analysis. 




C 

FIG. 2. Confidence contours of 68%, 95% and 99% in the (c, S7^)-plane for the case of marginal- 
izing h. The point in the figure, with the coordinate (0.21, 0.47), represents the best fit value, with 
Xmin = 173.44. Constraints from the age of the Universe give 0.95 < i^o^o < 1-05 (at the la 
confidence level), the dashed lines represent these two limits. 

Recent analyses of the age of old stars [35] indicate that the expansion time is in the 
range HGyr < t < 17Gyr at 95% confidence, with a central value t ~ 13Gyr. Following 
Krauss and Chaboyer [35] these numbers add 0.8 Gyr to the star ages, under the assumption 
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that star formation commenced no earlier than z — 6. A naive addition in quadrature 
to the uncertainty in H indicates that the dimensionless age parameter is in the range 
0.72 < H t < 1.17 at 95% confidence, with a central value H t ~ 0.89. More recent, 
Richer et al. [36] and Hansen et al. [37] found an age of 12.7 ± 0.7 Gyr at 95% confidence 
using the white dwarf cooling sequence method. For a full review of cosmic age see Ref . [4] . 
All in all, it seems reasonable to view ~ 12 Gyr to be a low limit of the cosmic age [38]. 
Now let us examine the age computation of the holographic dark energy model. The age of 
the Universe can be written as 

h=H « l L (i+%(z) ' (17) 

where Hq 1 represents the Hubble time with value Hq 1 = 9.778/i~ 1 Gyr. Using the best fit 
values, the holographic dark energy model gives the cosmic age t — 13.3 Gyr. This value is 
consistent with the above observational analyses. Furthermore, we impose a more rigorous 
test on it. The latest value of cosmic age appears in the Review of Particle Physics (PDG) 
[34] is given by a combined analysis of various observations [4], to = 13.7 ± 0.2. Using the 
data from Ref. [34], we can get the dimensionless age parameter range 0.96 ^ H to ^$ 1.05 at 
68% confidence, with a central value Hoto — 0.99. We also display the contours i^o — 0.96 
and H to = 1.05 (la) in Fig. 2. It can be seen explicitly that the la fit result of the SNe la 
data is almost excluded by this age test. 
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FIG. 3. Confidence contours of 68%, 95% and 99% in the (c, O^rJ-plane for the case of h = 0.64. 
The best fit values for the parameters are: Q? m = 0.41+g;^ and c = 0.46l°;25, with x^in = 175.90. 
Dashed lines represent the contours of cosmic age with to = 12, 13, 14, 15, 16 and 18 Gyr. 
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FIG. 4. Confidence contours of 68%, 95% and 99% in the (c, 0^)-plane for the case of h = 0.66. 
The best fit values for the parameters are: = °- 47 iaii and c = °- 21 ^aii' with Xmin = 173.44. 
Dashed lines represent the contours of cosmic age with to = 12, 13, 14, 15 and 16 Gyr. 
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FIG. 5. Confidence contours of 68%, 95% and 99% in the (c, 0^)-plane for the case of h = 0.71. 
The best fit values for the parameters are: = 0.49±0.03 and c = 0.05±0.02, with xL = 175.90. 
Dashed lines represent the contours of cosmic age with to = Hj 12, 13, 14 and 16 Gyr. 

Next, we will probe the sensitivity to the present Hubble parameter Hq in the analysis 
of the SNe data. We check that, keeping the parameter h fixed, how well the SNe and the 
cosmic age will constrain the holographic dark energy model. We fix h — 0.64, 0.66 and 0.71, 
respectively, and show the fit results in Figs. 3-5. From these figures, we notice that with the 
increase of the parameter h, c decreases evidently and the confidence contours in the (c, fi^J- 
plane shrink sharply. These figures clearly show that the SNe analysis is dependent on the 
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parameter h highly. Hence, finding out observational quantities which do not depend on H 
to jointly constrain the holographic dark energy model becomes very important and senseful. 
We also display age contours in these figures, but the age constraints are rather weak. The 
fits of SNe data indicate that the current SNe la data tend to support a holographic dark 
energy with c < 1, in other words, a quintom-type holographic dark energy. However, 
the authors of Ref. [17] tried to show another possibility. They re-examined the holographic 
dark energy model by considering the spatial curvature, and found that the holographic dark 
energy will not behave as phantom if the Universe is closed. However, the spatial flatness 
is a definitive prediction of the inflationary cosmology, and has been confirmed precisely by 
the WMAP. Thus we use the spatial flatness prior throughout this paper. 



The above analyses show that the supernovae data alone seem not sufficient to constrain 
the holographic dark energy model strictly. First, the confidence region of c — plane is 
rather large, especially for the parameter c. Moreover, the best fit value of is evidently 
larger. Second, it is very hard to understand the fit value of parameter c, 0.21, it seems 
odd because it leads to an unreasonable present equation of state, Wq = —2.64, the absolute 
value is too large. Third, even though the predicted cosmic age to — 13.3 Gyr is larger than 
the reasonable low limit of the cosmic age estimated by old stars, a more rigorous analysis 
implies that the fit result of the SNe la data is contradictive to the present data of cosmic 
age (dimensionless age parameter Hoto). Furthermore, our analysis shows that the fit of the 
SNe la data is very sensitive to the parameter Hq. Hence, it is very important to find other 
observational quantities irrelevant to H as complement to SNe la data. Fortunately, such 
suitable data can be found in the probes of CMB and LSS. 

In what follows we will perform a combined analysis of SNe la, CMB, and LSS on the 
constraints of the holographic dark energy model. We use a x 2 statistic 



where x|n i s given by equation (16) for SNe la statistics, Xcmb an d Xlss are contributions 
from CMB and LSS data, respectively. For the CMB, we use only the measurement of the 
CMB shift parameter [39], 



where z& cc = 1089 [3]. Note that this quantity is irrelevant to the parameter Hq such 
that provides robust constraint on the dark energy model. The results from CMB data 



IV. COMBINED ANALYSIS WITH CMB AND LSS 



z Z i Z I z 

X — XSN ' XCMB ' XhSS j 



(18) 




(19) 
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correspond to TZ = 1.716 ± 0.062 (given by WMAP, CBI, ACBAR) [4,40]. We include 
the CMB data in our analysis by adding Xcmb = [(^ ~~ 7^o)/W] 2 (see [41]), where 1Z is 
computed by the holographic dark energy model using equation (19). The only large scale 
structure information we use is the parameter A measured by SDSS [42], defined by 



where z\ = 0.35. Also, we find that this quantity is independent of Hq either, thus can 
provide another robust constraint on the model. The SDSS gives the measurement data 
[42] A = 0.469 ± 0.017. We also include the LSS constraints in our analysis by adding 
Xl S s = [(-A — A )/a A ] 2 (see [43]), where A is computed by the holographic dark energy 
model using equation (20). 

Note that we have chosen to use only the most conservative and robust information, 1Z 
and A, from CMB and LSS observations. These measurements we use do not depend on 
the Hubble parameter Hq. Furthermore, by limiting the amount of information that we use 
from CMB and LSS observations to complement the SNe la data, we minimize the effect 
of the systematics inherent in the CMB and LSS data on our results. Figure 6 shows our 
main results, the contours of la, 2a, and 3a confidence levels in the c — plane. The 
best fit values for the model parameters are: h = 0.65, fi^ = 0.28 ± 0.03, and c = 0.8llo;iJ, 
with Xmm = 176.67. We see clearly that a great progress has been made when we perform 
a joint analysis of SNe la, CMB, and LSS data. Note that the best fit value of c is also less 
than 1, while in la range it can slightly larger than 1. Now the fit value of is roughly 
as the same as that of the ACDM model (WMAP result), but h is still slightly smaller. It 
should be pointed out that a slightly lower value of h is, however, in agreement with the 
observations of [44] which can accommodate lower values of h ~ 0.6. We also see from the 
figure that the fit result of the joint analysis is consistent with the dimensionless age range 




(20) 



0.96 < H t Q < 1.05. 
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FIG. 6. Results for analysis of SNe+CMB+LSS data, la, 2a, 3a confidence levels in the c— Sl^ 
space, marginalizing over h. The best fit values are: h = 0.65, n° m = 0.28, and c = 0.81, with 
Xmm = 176.67. Cosmic age test, the dimensionless age parameter range 0.96 < H^t® ~ 1.05, also 
shown. 

We now compare the fit results of the SNe analysis and the joint analysis of SNe, CMB, 
and LSS, and discuss the different cosmological consequences. Let us first look at the 
likelihood distributions of the parameter c in the two fits. In Fig. 7 we plot the 1-dimensional 
likelihood function for c, marginalizing over the other parameters. The big difference can lead 
to rather different conclusions for some important cosmological parameters, such as today's 
equation-of-state parameter of dark energy, today's deceleration parameter, etc. For the 
best fit results of the two fits, we plot the evolution behaviors of the deceleration parameter 
q and the equation-of-state parameter of dark energy w in Fig. 8. For the SNe+CMB+LSS 
joint analysis, from the figure, we see that the deceleration parameter q has a value of 
go = —0.61 at present. The transition from deceleration to acceleration {q{zx) = 0) occurs 
at a redshift of zt = 0.63. The equation-of-state parameter w is slightly smaller than —1 at 
present, wq = —1.03. For comparison, we list these values for the alone analysis of SNe data: 
qo = —1.60, zt = 0.27, and wq = —2.64. Therefore, obviously, the results of the combined 
fit seem more reasonable. From a joint analysis of SNe+CMB+LSS data, one may obtain 
within the framework of the holographic dark energy model a fairly good idea of when the 
Universe began to accelerate and how fast the present acceleration is. By contraries, the 
cosmological consequences given by the alone analysis of the SNe data seem unreasonable. 
Comparing our plots in Fig. 8 with the model-independent plots in Ref. [27] (which also use 
data from [33]), we find that the holographic plot for c = 0.81 case is in good agreement 
with those model-independent plots for the redshift range z = — 2, while the c = 0.21 case 
dose not accord. Moreover, it should be mentioned again that these results demonstrate 
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that the best fit to SNe+CMB+LSS observations also favors a quintom-type holographic 
dark energy. 
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FIG. 7. Likelihood distributions of parameter c in the fits of SNe only and SNe+CMB+LSS. 
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FIG. 8. Deceleration parameter q and equation of state of dark energy w, versus red-shift 
z, from two fits, SNe only and SNe+CMB+LSS. Solid lines correspond to the joint analysis of 
SNe+CMB+LSS, with the parameters c = 0.81, = 0.28, and h = 0.65. Dashed lines correspond 
to the alone analysis of SNe, with the parameters c = 0.21, = 0.47, and h = 0.66. 
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V. EXPECTED SNAP ANALYSIS 



Finally we consider how well the proposed Supernova/Acceleration Probe (SNAP), may 
constrain the parameters c and Q^. The SNAP mission is expected to observe about 2000 
type la SNe each year, over a period of three years, according to the SNAP specifications. To 
find the expected precision of the SNAP, one must assume a fiducial model, and then sim- 
ulate the experiment assuming it as a reference model. We will use SNAP specifications to 
construct mock SNe catalogues [45]. Following previous investigations [46], we assume, in 
our Monte Carlo simulations, that a total of 2000 supernovae (roughly one year of SNAP 
observations) will be observed with the following redshift distribution. We consider, 1920 
SNe la, distributed in 24 bins, from z = to z = 1.2. From redshift z — 1.2 to z — 1.5, 
we assume that 60 SNe la will be observed and we divide them in 6 bins. From z — 1.5 
to z — 1.7 we consider 4 bins with 5 SNe la in each bin. All the supernovae are assumed 
to be uniformly distributed with Az = 0.05. To fully determine the x 2 functions, the error 
estimates for SNAP must be defined. Following [45], we assume that the systematic errors 
for the apparent magnitude, m, are given by 



which are measured in magnitudes such that at z — 1.5 the systematic error is 0.02 mag, 
while the statistical errors for m are estimated to be a sta = 0.15 mag. We add both kinds 
of errors quadratically 



where n» is the number of supernovae in the i'th redshift bin with width Az 0.05. Now 
let us assume a spatial flat ACDM model as a fiducial model, and analyze the holographic 
dark energy model fit. We aim to show if the Universe is indeed described by the ACDM 
model, how well the the fitting of a holographic dark energy model to be. 
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FIG. 9. Predicted confidence contours of 68%, 95% and 99% in the (c, fi^)-plane for the SNAP 
mission are shown. We consider a fiducial model with $7^ = 0.27 and h = 0.71. The best fit values 
for the parameters are: fij^ = 0.23 and c = 0.92. 

In Fig. 9 we display the results of our simulation assuming a ACDM model as fiducial 
model with = 0.27 and h = 0.71. In the fit, we marginalize over the nuisance parameter 
h. The best fit values for the model parameters are: = 0.23 and c = 0.92. From this 
figure it is clear that SNAP mission is able to place rigorous constraints on the holographic 
dark energy model. On the other hand, we notice with interest that, even though the dark 
energy in the Universe is exactly described by a cosmological constant A, the precision type 
la supernova observations will still support a quintom-type holographic dark energy. 

VI. CONCLUDING REMARKS 

In this paper, we investigated constraints on the holographic dark energy model from 
current and future SN la observations. We considered a spatially flat FRW Universe with 
matter component and holographic dark energy component. For the holographic dark energy 
model, the numerical constant c plays a very important role in determining the evolutionary 
behavior of the space-time as well as the ultimate fate of the Universe. In a holographic 
dark energy dominated Universe, the case of c = 1 corresponds to an asymptotic de Sitter 
Universe; The choice of c < 1 will lead to dark energy behaving as quintom, and in this case, 
the Gibbons-Hawking entropy will eventually decrease as the event horizon shrink such that 
violates the second law of thermodynamics; While the case of c > 1 does not violate the 
second law of thermodynamics, in this situation, the evolution of the corresponding space- 
time avoids entering de Sitter phase and Big Rip phase. Though the choice of c = 1 is favored 
theoretically, other possibilities can not be ruled out from the viewpoint of phenomenology, 
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and only experiments and observations are capable of determining which choice is realistic. 
We derived model parameter ranges from the analysis of the present available SNe la data, 
and then imposed a rigorous test using the new data of age of the Universe on the derived 
parameter region. However, only the supernova analysis seems not sufficient to be able to 
precisely determine the value of c. For improving the result of the analysis, we perform a 
joint analysis of SNe, CMB, and LSS data to the holographic dark energy model. 

The results of the SNe analysis show that the holographic dark energy model behaves as 
quintom in la confidence level, consistent with the current SNe la data. However, when we 
perform a rigorous age test on this analysis result using the latest PDG data, the la allowed 
range is almost ruled out. Moreover, in this case, the confidence regions in the parameter- 
plane are rather large. We also probe the sensitivity to the present Hubble parameter H 
in this analysis. We do this by fixing the parameter h, and find the allowed regions of the 
parameters shrink sharply as h increases. Hence, a combined analysis of SNe data with other 
observational data becomes important. A joint analysis of SNe+CMB+LSS produces more 
reasonable results: c = 0.81, = 0.28, and h = 0.65, leading to the present equation of 
state of dark energy is w = —1.03, and the epoch at which the Universe began to accelerate 
is zt = 0.63. The confidence regions in this analysis become more compact. On the whole, 
the analysis indicates that the case of c < 1 is consistent with the present SNe la data. 

We are also interested in the constraints on the holographic dark energy model from a 
fictitious future supernova experiment. An expected SNAP fit, by using the ACDM model 
as fiducial model to generate mock observational data, shows that the case c < 1 (c = 0.92) 
also favored. Obviously, a large number of supernovae at high redshifts, as well as better 
knowledge of the values of H and are therefore required to draw firm conclusions about 
the property of the holographic dark energy. We expect that a more sophisticated combined 
analysis of various observations will be capable of determining the value of c exactly and 
thus revealing the property of the holographic dark energy. 
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